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SUMMARY 

Sources of high-energy electrons are encountered in space (the magnetically trapped Van Allen belt 
electrons) and in ground-level high-energy physics laboratories. It is important to be able to predict 
the damage resulting to human beings and radiation-sensitive equipment near these sources. In this 
report, various techniques for the calculation of electron and bremsstrahlung dose deposition are des- 
cribed. New energy deposition, transmission, and reflection coefficients for electrons incident on plane 
slabs for angles of 0, 30, 60, 75, and 89.9 deg and energies of 0.5, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, and 
10.0 MeV are presented, and methods for their use in electron dose calculations are developed. A method 
for electron dose calculations using the M straight-ahead ,, approximation is also developed, and the 
various methods are compared and found to be in good agreement. Accurate and approximate methods of 
calculating bremsstrahlung dose are derived and compared. The approximation is found to give good order 
of magnitude estimate of dose where the electron spectrum falls off exponentially with energy. The pri- 
mary weakness of both calculations is the approximation of the bremsstrahlung source angular distribution* 
the actual distribution is not easily determined. More work needs to be done in this area 


USE OF BERGER'S MONTE CARLO DATA IN 
ELECTRON DOSE CALCULATIONS 

Probably the most successful attack on the 
electron transport problem has been by Martin 
Berger U~3] using Monte Carlo methods and multiple 
scattering techniques. His present set of pro- 
grams will take an incident beam of electrons or 
photons and follow both the primaries and any 
secondary electrons or photons produced. The out- 
put includes almost any quantity of interest de- 
pending on the program option used. Because the 
simulation is so thorough, the program is very 
complex and requires large amounts of computer 
time on one of the larger machines available to 
complete an accurate calculation. Thus, the pro- 
gram's primary usefulness is in generating basic 
data for incorporation into other programs using 
more simplified approaches to the problem. 

Berger's program was made available to 
Marshall Space Flight Center (MSFC) , and a number 
of calculations have been performed here for com- 
parison with experimental workf^], As a second 
study systematic calculations covering the energy 
range of from 0.5 to 10.0 MeV for a number of 
angles of incidence by a monodirectional beam on 
a plane aluminum slab were made. The geometry is 
shown in Fig. 1. Of particular interest were the 
electron energy and number transmission and re- 
flection and the energy deposition coefficients, 
because they can be incorporated into an electron 
dose or number deposition calculation involving an 
arbitrary incident electron energy and angular 
distribution. To minimize the computer time used 
no photons were followed, and only enough electron 
histories (2500 through 7000, depending on the 
angle of incidence) were sampled to get good sta- 
tistics for the reflection, transmission, and 
deposition coefficients. (A much higher number of 
histories would have been required to get good 
statistics on one of the differential quantities 
such as exiting energy spectrum.) The fact that 
secondary photons were not followed means that 
there is a slight underestimate both in the number 
of coefficients due to missing tertiary electrons 
produced by the photons and in the energy factors 
due to energy transported by secondary photons. 
Because over most of the energy range of interest 
electrons lose only a small fraction of their 
energy as photons, both these effects should be 
small. The radiative yield, the fraction of an 


electron's energy lost as bremsstrahlung in stop- 
ping, is 0.003324 at 0.25 MeV and 0.07721 at 10.0 
MeV in aluminum^ . This quantity should set an 
upper bound on the possible error induced by the 
limitations imposed at least for the energy factors. 



FIGURE 1: BASIC GEOMETRY USED FOR CALCULATIONS 
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There are several approaches for attacking 
a dose deposition calculation using the results 
tabulated by Berger’s electron transport program. 
The most obvious is to use the internal energy 
spectrum and instantaneous stopping power. Two 
less obvious but more efficient methods - one 
reasonably exact and the other approximate - will 
be described here. 

The exact calculation makes use of a quantity 
Berger calls energy deposition (Figs. 2 and 3), 
which is the average energy deposited per unit mass 
per electron at a given depth into a plane, in- 
finitely thick slab by, in this case, a monodirec- 
tional beam of electrons. [This does not follow 
the initial model (Fig. 1) because there is 
material behind the receiver point.] Berger calcu- 
lates the energy deposition by breaking up the slab 
into thin layers and keeping an inventory of energy 
deposited in each layer and then dividing this 
quantity by the incident number of initial electron 
histories and layer thickness. Thus, 

N. 


P(E, 0, X.) = 

where p is the energy deposition in MeV/g/unit 
current, E is the incident energy in MeV, Q is 
the incident beam angle measured from the normal to 
the slab, Xj is the depth to the center of the jth 
layer, Nj is the number of electrons penetrating 
layer j , AE^ j is the energy in MeV deposited by 
the ith penetrating electron in the jth layer, N 0 
is the number of initial electron histories, and 
AXj is the thickness of the jth layer in g/cm 2 . In 
Berger’s tabulation Xj is measured in fractions of 
an electron pathlength at the incident energy. The 
electron mean pathlength is the average length of 
the zig-zag path followed by an electron in 
stopping as opposed to the mean range which is the 
average straight-line distance traversed. It is 
given by: 

■ < 2 > 

where r Q (E) is the pathlength at energy E and 
S t (E’) is the total instantaneous stopping power. 
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If there is an angular and energy flux distri- 
bution given by <I> 0 [E,S(0, <j>) ] with electrons 
incident per unit at E and per unit solid angle in 
the direction Q, then the incident current as used 
by Berger is 


j 0 (if) = * 0 (E, n) cos (e) , o) 

and the dose at a depth Z in g/cm^ is given by: 


D( Z) = K J f pi E,0,Z/r o (E)] * 0 <E,fi) 

cos 0 dE dlf 


7? E 


(4) 


where K is a units conversion constant. Dose is 
more often measured in rads than MeV/g, in which 
case K = 1.60 x 10"® rads/ (MeV/g) . * In the case 
where there is a monodirectional beam incident at 
an angle 0 from the slab normal, Equation (4) 
becomes : 

D( Z, 0) = K cos 9 f p[E,6, Z/r 0 (E)] (5) 

E 

4> 0 (E) dE . 


Another case of interest is that in which the 
distribution is half-space isotropic. In this 

case the dose is given as 
7 r 

2 

P i 80 (Z)=K f f PlE.0, Z/r 0 (E)] (6) 

E 0 

cos 0 sin 0 do $o (E) dE . 

Thus, a half-space energy deposition function may 
be defined by 

T 

P iS0 (E,X)= / P(E,0,X) (7) 

cos 0 sin 0 d 9 

*To convert rads to SI Units in joules per 
kilogram, multiply rads by 0.01. 
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This function is tabulated in Table 1. If Equation 
(7) is used, Equation (6) becomes: 


How these factors may be used in energy depo- 
sition or dose calculations will now be considered. 
From conservation of energy. 


d . 8 o ( Z)=K f P lBQ [E, Z/r,(E)J ( 8) 

E 

*0 (E) dE . 


E =■ E [Ag (E.O.X) + U E (E.9.X) 

+ T e (E.O.X)] (10) 


The approximate method for calculating elec- 
tron dose makes use of the energy transmission and 
reflection factor of Berger to derive an approxi- 
mation to the energy deposition function. The 
electron energy current transmission and reflection 
and number current transmission and reflection 
factor shown in Figs. 4 through 11 are defined as 
follows : 


or 

1 = A e (E, 0.X) + U E (E.O.X) 
+ T e (E.O.X) 


( 11 ) 


T n (E.O.X) 


N t (E.O.X) 
Jo (0) 


A n (E.O.X) = 
T e (E.O.X) = 


Ae (E.O.X) = 


N r (E.O.X) 

J*(0) 

N_ 

T 

Z E Ti (E.e.X) 

;=l 

E J 0 (0) 

n r 

E E r . (E.O.X) 
E J 0 (0) 


are obtained, where Ue(E, 0,X) is the fraction of 
the energy either deposited in the slab or radiated 
as bremsstrahlung. Since in the energy range of 
interest the radiated component is small, it will 
be assumed negligible. If a quantity f(X T ,E,0,X) 
dX f is defined, which is the fraction of the 
incident energy deposited between X f and X 1 + dX* , 
then 

X 

U E (E.O.X) = / f (X\ E.O.X) dx' . ( 12 ) 

( 9 ) 0 

Using Equations (11) and (12) yields 

X 

/ f(X* E.O.X) dX’ = 1 -A_ (E.O.X) 

0 ** 

(13) 

-T E (E.O.X) . 


where T^ is the number transmission factor; Nj is 
the number of electrons passing through a slab X 
thick; An is the number reflection factor where 
Nr is the number of electrons reflected from a 
slab X thick; Te is the energy transmission factor 
where E?i is the energy of the ith transmitted 
electron; and Ae is the energy reflection factor 
where Eri is the energy of the ith reflected 
electron. 
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K I CURE 4: ENERGY TRANSMISSION COKrT 1 Cl ENT AT 0.5, J.O, AND 
10. 0 MeV .-OR A NORMALLY INCIDENT ELECTRON UK AM 


Taking the derivative with respect to X gives 

f(X. E.0.X) + J -g- (X', E.0.X) dX* = 

0 

d (14) 

"dX [A E (E ’ 0 ’ X) + T E (E,0 t X] . 



ELECTRONS INCIDENT AT VARIOUS ANGLES 
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ENERGY REFLECTION COEFFICIENT A E ENERGY REFLECTION COEFFICIENT A £ 







Z/r 0 (E) 

FIGURE 10: NUMBER REFLECTION COEFFICIENT AT 0.5, 

2.0, AND 10.0 MeV FOR A NORMALLY INCIDENT ELECTRON 
BEAM 



ELECTRONS INCIDENT AT VARIOUS ANGLES 


The approximation in this method involves ignoring 
the integral term on the left in this equation. 
Since f(X',E,0,X) increases with increasing X (due 
to reflected electrons), the derivative is always 
greater than or equal to zero. Thus the integral 
is positive, and ignoring it gives a conservative 
estimate of f(X,E,0,X). It is difficult to justify 
this approximation except to observe that in 
practice it yields results comparable with those 
of the previously described method. It was deve- 
loped because Berger's older publication presented 
only the transmission and reflection factors, the 
energy deposition factor having become available 
only recently. 

The energy deposition function (in units of 
MeV/g) is given by: 


p(X,E,0,X) 


Ef(X t E,e,X) 
r o (E) 


(15) 


-E d 


P(X ’ E ’*’ X) ~ rTTffdx tV E ’*> x ) 

+ T e (E, 0 ,X)] (16) 


It is interesting to note that f(X',E,0,X) for 
X' < X always has a contribution caused by re- 
flection from portions of the shield beyond X* but 
that f(X,E,0,X) does not. Thus it approximates the 
quantity needed for an energy deposition cal- 
culation in our original geometry. 

Because of the definition of p in the deriva- 
tions the material of the receiver must be the 
same as that of the shield. An approximate 
correction for estimating the dose for a different 
receiver can be made by multiplying the single 
material calculation by the ratio of the collision 
stopping power in the shield material at some 
typical energy for the exiting electron spectrum. 
Fortunately, the energy selected is not particul- 
arly important since the ratio of two electron 
stopping powers is not a sensitive function of 
energy except at very low energies ( > 0.01 MeV). 
Some estimate of the accuracy of the approximation 
can be found by observing the variation of the 
ratio with energy. A particularly interesting 
case is that of an aluminum shield and a tissue 
receiver. This case, using Berger's stopping 
power datat3] is given in Table 2. 

Taking the simple average of ratios between 
energies of 0.1 to 10.0 MeV yields a correction 
(1.30) that will be within about 3% of an exact 
calculation under most conceivable conditions. 

(For typical incident energy spectra encountered, 
it would be unusual for the average energy of the 
exiting spectrum to be less than 0.04 MeV). 

Thus, for an aluminum shield and a tissue receiver 
Equations (5) and (6) become, respectively: 

D(Z,0) = 1.3 K cos e f p [E, 0, Z/r 0 (E) ] 

E 

*0 (E) dE 


and 

D. so ( Z) = 1.3 K f p. so |E,Z/r 0 (E)] (18) 

E 
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TABLE 2: RATIO OF TISSUE COLLISION STOPPING POWER TO ALUMINUM 

COLLISION STOPPING POWER AS A FUNCTION OF ENERGY 


Energy 

(MeV) 

_ dE/dX tissue 
dE/dX aluminum 

------ X 100 

(%) 

0.01 

1.383 

6.0 

0.02 

1.350 

3.7 

0.04 

1.325 

1.9 

0.06 

1.313 

1.0 

0.08 

1.305 

0.4 

0.1 

1.300 

0.0 

0.2 

1.285 

-1.2 

0.4 

1.276 

-1.9 

0.6 

1.272 

-2.2 

0.8 

1.271 

-2.3 

1.0 

1.272 

-2.2 

2.0 

1.281 

-1.5 

4.0 

1.301 

0.1 

6.0 

1.317 

1.3 

8.0 

1.329 

2.2 

10.0 

1.340 

3.0 


Curve fits of P(E,0,X) have been found to be 
very useful for computational purposes. Fit over 
X of the form M 


P(E, 0,X) = e 



(19) 


has been found to give satisfactory results for 
normal incidence and half-space isotropic inci- 
dence, two cases of special interest. Fits of 
the above form were made and then the coefficients 
were fit as a function of energy. For normal 
incidence the coefficients are given by: 


A - non °* 96 3E 

- 0* 913 e + 0.021E + 0.215 

A 2 = 5. 0 - 0. 491E 


a 3 

*4 


57. 573 (E - 5. 0)/(E + 29. 98) 
-i.6E°- 837 


( 20 ) 


A 1 


0.52 + 0. 09854 


e -1.4U8 


A 2 

A S 

A 4 


.-°' 821E - 1.0 

3.253 e -°- 323E *5.8 


( 21 ) 


A 5 - -15. 4375 + i.;5542E - 0. 0786077E 2 


Because the shape of the energy deposition 
function for angles near 90 deg is different from 
that near 0 deg, a good fit by a single functional 
form is difficult to achieve over the whole range 
of directions. Best results will probably be 
obtained by interpolation from the actual data or a 
smoothed set derived from it. The function f(E,0, 
X) derived using the reflection and transmission 
factors was also fit but by a different form. The 
fit was to the form 


and for half-space isotropic the coefficients are 
given by: 


1 -T e <E,0,X) - A e (E,0,X) = 
A [ t _ e (BX+ CX 2 + DX 3 )] 


( 22 ) 
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where A, B, C, and D are functions of E and 0. 

Then, taking the derivative with respect to X and 
using Equation (14) yields 

f(E,0,X) = (23) 

, (BX+ CX 2 + DX 3 ) 

-A (B + 2CX + 3DX Z ) e 

Fits have been derived for the two cases of special 
interest. For a normally incident beam (0=0 deg) 
in the energy range 0.5 to 10.0 MeV, 

0.893 + 1.682E 
A ' 1.0+ 1.665E 


11. 0E 

B " 11. 0E +6.0 
C = 4.2 exp(-0.47E) 

and 

D = 5. 16 . 

Berger has published a transmission and reflection 
coefficient for a cosine law source (half-space 
isotropic flux)t^l. In this case fi so (E,X) is a 
relatively insensitive function of E, and satis- 
factory results are obtained by taking an average 
curve for all energies. For energy deposition. 


A = 0.439 


B - -2.08 


C = -3.54 

(25) 

D = -6.08 


The electron pathlength r 0 (E) in aluminum taken 
from Berger [2] has also been fit by the following 
form: 

r 0 (E) = 

(26) 

(1.33 - 0.019E) 

(■s/o.2713E 2 + 0.0121 -0.il) , 


which is within 2% of actual curves for energies 
greater than 0.3 MeV and within 5% between 0.2 and 
0.3 MeV. 


ELECTRON TRANSPORT USING THE STRAIGHT-AHEAD 
AND CONTINUOUS SLOWING DOWN APPROXIMATIONS 

A method commonly used for describing the 
transport of protons is to assume that the 
particle travels through the shield along its 
incident direction losing energy continuously 
according to some stopping power law. Thus, 
according to the approximations, the energy and 
direction of the particle at any point in the 
shield is completely predictable. This method 
has the advantage of providing an energy spectrum 
at the internal point of interest that may be 
used to determine such things as secondary pro- 
duction sources as well as to calculate dose 
deposition. Its disadvantages for application 
to electron transport are twofold: First, 
electrons are more likely to be scattered from 
their original direction than protons, and 

*This mean range is not Berger’s mean pathlength. 


second, electrons can suffer large energy losses 
in a single interaction. Thus the straight-ahead 
approximations should not be expected to be es- 
pecially applicable. The approximations have been 
applied with some success, however. 

Since the validity of the approximations in 
applications to electron transport are somewhat 
questionable, one would prefer to be on the con- 
servative side in any estimate of particle energy 
or number. By using the extrapolated electron 
range rather than the mean range* for definition 
of other required quantities this can be accom- 
plished. The extrapolated range is defined as 

shown in Fig. 12, where it can be seen that few 
electrons penetrate beyond this distance. A very 
good fit to the extrapolated range in aluminum is 
given by: 

r ex =H< E) =y^- b (27) 

over the energy range from 0.0 to 16.0 MeV. For 
R(E) in g/cm^ E is in MeV, a is 1.92, and b is 
0.1115], 



FIGURE 12: DEFINITIONS OF ELECTRON MEAN RANGE r M£AN » 

EXTRAPOLATED RANGE r^, AND MAXIMUM RANGE 


Given the electron range, one can determine 
the relationship between the initial energy and the 
energy after passing through thickness Z of 
material as follows : 

Since 

R (E) = R (E') + Z , (28) 

where E is the initial energy and E' is the energy 
at depth Z, 
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E = g(E',Z) =R _1 [R(E') + Z)]; (29) 


R _1 denotes the inverse of the function R. Thus, 
using Equations (27) and (29) gives 


E = g(E\z) = 



The relationship between the external differential 
flux $> q (E) and the flux at depth Z, <J> Z (E'), is 
given by: 


Equation (36) is derived for normal incidence. For 
a beam incident at angle 0 to the normal, Z is 
replaced in the equation by the slant distance to 
the dose point, Z/cos 0, and for half-space iso- 
tropic flux the dose is given by: 

7T 

2 

D. (Z)=K f f [g< E, » Z/ C ° s 0 )] 
iso i 

(37) 


COMPARISON OF THE THREE METHODS FOR 
DOSE CALCULATION 


^ J+ b 2 + Z/cos e 
g(E\Z/pos 0 ) 


dE* sin 0 d 0 


* z (E 1 ) = [g(E\Z)] 


dE’ 


(31) 


The derivative enters the equation because of 
change of* energy e to E* units; or as one can see, 
lower energy electrons lose energy faster, thus 
changing particle densities. Using Equation (30) 
gives 


* J T-l 1 


‘dE 




+ b 2 + Z 




"12 


+ b z + Z 


( 32 ) 


The electron dose deposited at depth Z is 
given by: 


D(Z) =K / * (E') S(E') dE' (33) 

E 


Typical electron spectra encountered in 
practical applications are exponential in nature. 
Figures 13 through 18 show comparisons of the 
three methods for spectra of the form 

$0 (E) = Pe PE ( 38 ) 

where P varies from 0.25 to 6.0. The energy inte- 
gration limits are 0.0 to 20.0 MeV. The agree- 
ment among the three methods is fairly good con- 
sidering the magnitude of the attenuation that 
occurs. The two methods using Berger 1 s data have 
significant disagreement only for high shield 
thicknesses and isotropic spectra. The disagree- 
ment here is probably because a significant frac- 
tion of the dose is from electrons with energies 
above 10 MeV where both methods are using extra- 
polation on the data. The disagreement between 
the methods using Berger 1 s data and the straight- 
ahead method are somewhat larger, but the errors 
still are not so large as to cause question about 
the validity of any one of the three approaches. 


where S(E') is the instantaneous collision stopping 
power in the receiver. One can derive an approxi- 
mation for the stopping power by using the deriva- 
tive of the range 


S(E’) = - 


dR(E') 

dE' 


««•»-# m*' 


(34) 


(35) 


By a fortunate accident the approximate stopping 
power derived from aluminum extrapolated range 
data is a good fit to tissue collision stopping 
power. (Tissue is the receiver most commonly 
used.) The fit is within 5% of Berger's tabulated 
dataL2] i n the interval 0.15 to 4.5 MeV and with- 
in 13% from 0.08 to 10.0 MeV. 

Combining the results of Equations (31), (32), 
(33), and (35) and simplifying the tissue dose 
behind an aluminum shield Z thick yields 
D(Z) = 


K / * 0 [g<£\Z)] 

n * 




g(E',Z) 


(36) 

dE'. 



COMPARISON OF ELECTRON TISSUE DOSES CALCULATED 
USING BERGER’S ENERGY DEPOSITION DATA AND 
TRANSMISSION - REFLECTION DATA FOR A NORMALLY 
INCIDENT SPECTRUM OF THE FORM i (t) . p eXD 
{-PE)e/cn* - MeV ° P 
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section differential in photon energy in units 
of photons/g-MeV-(e/cm2) and (E*S) is the 
electron energy and angular distribution at depth 
Z in units of e/cra 2 -MeV-sr. 

The macroscopic cross section is given by: 

Z'V-l'-Sr • <«> 

y 

23 

where N 0 is Avagadro’s number 6.02 x 10 atom/ 
mole, A is the atomic weight of the material at 
the source point in g/mole, and da/dE is the 
microscopic bremsst rah lung production cross 
section differential in photon energy in units of 
photons /atom . (e/cm 2 ) . MeV. The microscopic 
cross section presently used for the MSFC cal- 
culation is given by: 



tttMPAKlSJK or U.KCTKOX TISSLL DOSI.S WLULXTI D 

IS I XI. m; ri;i:r*s tiuxs.mission-rkki.kltiox data axd mr 

STRAIGHT- AH HAD MOUKL KUK A HAI.K-KPAIK 1SHTKDIMC 
IXCIDKXT SI’HIJ I'Rl'N OK THL FORM ; (I ) . )• 

Wcra- - MoV 11 1 


BREMSSTRAHLUNG DOSE CALCULATIONS 

Previously, consideration has been confined 
to the dose deposited by electrons; this is satis- 
factory for thin shields. However, when the 
shields are thick enough to remove a large fraction 
of the primary electrons, dose deposition by 
bremss trahlung must be taken into account. (The 
point where this occurs is usually less than 3.0 
g/cm2.) To do a bremss trahlung dose calculation 
one must generate a bremss trahlung source distri- 
bution and then transport the bremss trahlung from 
this source to the dose point. To generate the 
bremss trahlung source distribution, one needs the 
electron energy and angular distribution at the 
source point and the bremss trahlung production 
cross section. 

Because the interior electron angular distri- 
bution is not easily determined and since the 
bremss trahlung production cross section differen- 
tials in angle are not particularly accurate, a 
simplification commonly made is to use cross 
section differential in photon energy only and to 
make some assumption about the bremss trahlung 
source angular distribution. (This assumption 
will be examined later.) In this case, the 
bremss trahlung differential energy flux source or 
depth Z' and photon energy Ey is given by: 

S(E,z') = 

y 

/ / E Tj .E) V (E.n) dEdjf 
E ~ 7 r L (39) 

in units of MeV/(MeV-g), where £ (Ey,E) is the 
macroscopic bremsstrahlung production cross 


where 



(41) 



The terms Hq and H are the initial and final total 
electron energies in me 2 units given by 



E - E 


H = 


me 




Pq and P are the initial and final electron 
momenta given by 

P 0 = >7 Ho 2 - 1 


= s/h 2 - 


and K is the photon energy or momentum in me 2 or 
me units, respectively. 


E 

K = 

me 


811 




(47) 


The term Z Q is the atomic number of the source 
material, r Q is the classical electron radius 
2,82 x 10“!3 cm, me 2 is the rest mass energy of 
an electron, and F E (E,Ey) and C(E) are two cor- 
rection factors. Except for the correction 
factors this is formula 3BN from Koch and Motz , 
who give a complete description of the cross 
section. The correction factor F E (E,E^) given by 


F E (E ’V 


where 




E < 2. 0 MeV 

5t> 0.01 <«> 

E 


1.0 


otherwise , 



and 


0 = 



me* 


- E + me 

y 


(43) 


(44) 


is the Elwert nonrelativistic coulomb correction, 
and the correction C(E) is an empirical screening 
correction as shown by Koch and Motzl^l. Koch and 
Motz further discuss both these corrections l 6 J . 

The prime advantage of the straight-ahead 
approximation is that it does provide an electron 
spectrum at a given depth. Thus, using the result 
derived, for a normally incident electron beam, 
the bremsstrahlung source is given as 


w z >- 


/ E £(E ,E> [g(E, Z*)] 

E y 7 


dgtE.Z 1 ) 

dE 


(45) 
dE , 


and for a half-space isotropic incident electron 
distribution the bremsstrahlung source is 



(46) 


* o [g(E,z'/cos0)] dS ^ E ’^ C ° S ^ ^ sin 0 d0 , 


where 4> 0 , g, and dg/dE are as defined in the third 
section. 

The bremsstrahlung transport and dose cal- 
culation is relatively straightforward for a 
plane slab geometry given the source energy dis- 
tribution and some source angular distribution. 
Using point kernel attenuation with dose build- 
up factors, the dose is given by: 


Z 

D( Z) = K f f M E 
E 0 

y 


<E ) S(E,z') B(Ey, Z, Z*) 


Gfri (E )|Z-z'|] dE dZ f , 

my y 


where y E (Ey) is the gamma ray energy absorption 
coefficient for the receiver material; Z is 

the shield thickness; B(E,Z,Z f ) is a dose buildup 
factor, depending on the source angular distribu- 
tion chosen; and G[y m (E T ) |Z-Z T | ] is the attenuation 
kernel, also depending on the source angular dis- 
tribution chosen where y m (E y ) is the gamma ray mass 
absorption coefficient. 

Much .o# the wide disagreement among various 
bremsstrahlung calculations can be traced to 
assumptions made about either the incident electron 
angular distribution expected in the problem or 
the angular distribution of the bremsstrahlung 
source. Perhaps the least conservative assumption 
about the source distribution that can logically be 
made for deep penetrations is that it is isotropic, 
and the most conservative assumption is that all 
the photons are emitted normally into the slab. In 
the first case the attenuation kernel is given by: 


Gbi m ( V |Z - Z, U = 


E i ^ m (E y )lz " z,|] 
2 


(48) 


where El is the first exponential integral and 
plane isotropic buildup factors are usedL»J . In 
the second case the attenuation kernel is 


-M (E ) |Z-Z f | 

G ^ m <V IZ - Z ' ll = e m T • (49) 

and plane monodirectional buildup factors are 
used. Goldsteint 7 ] tabulates these for infinite 
media, which should give a conservative estimate 
of the dose. These should be used cautiously, 
however, because they do not extend low enough in 
energy and extrapolation is dangerous. Figure 20 
gives a comparison of these two cases for a half- 
space isotropic electron spectrum of exponential 
form Incident on an aluminum shield with a water 
receiver (simulating tissue). Plane Isotropic 
buildup factors were used in both calculations so 
that the difference observed is caused by the 
attenuation kernel. Actually, if correct buildup 
factors were used, the normal incident source 
case would be slightly lower. 


Because of the built-in bias of using the 
extrapolated range in the straight-ahead approxi- 
mation for calculating the source distribution, it 
is felt that there is no need to use the most con- 
servative source angular distribution. Instead it 
is assumed that the bremsstrahlung source is half- 
space isotropic toward the receiver by the atten- 
uation function given in Equation (48) multiplied 
by a factor of two. Still, the plane isotropic 
buildup factors W are used, making the calculation 
slightly more conservative. Figure 20 shows the 
results of an MSFC bremsstrahlung calculation for 
a half-space isotropic electron spectrum of expon- 
ential form incident on an aluminum shield with a 
water receiver (simulating tissue) . Comparing 
Figs. 19 and 20 reveals that the MSFC calculation 
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(Fig. 20) is the most conservative of the three up 
to 4 or 5 g/cm 2 . Above that depth, the mono- 
directional source calculation is more conservative 
by as much as a factor of three compared to the 
MSFC results . 




M^lKI 20: BHKMSSTKAHLUNt DUSE CALCULATIONS ASSUMING A MALE-SPACE 

" ISOTROPIC SOURCE FOR A HALF-SPACE ISOTROPIC INCIDENT 

ELECTRON SPECTRIN OF THE FORM ; o (E) - |> uxp<-PK)Wcm- - 'k-\ 


Since the bremsstrahlung dose calculational 
method described so far involved triple or quad- 
ruple numerical integration and since the functions 
integrated are exceptionally difficult to inte- 
grate, the calculation can be performed convenien- 
tly only on a relatively large computer. In many 
cases, all that is needed is a rough order of 
magnitude estimate useful only in determining if 
a problem exists. For space applications, where 
most of the electron spectra encountered fall off 
exponentially with increasing energy, a conserva- 
tive estimate of bremsstrahlung dose can be cal- 
culated as follows : 

First, assume that all the electrons penetrate 
to the source plane and that at that plane their 
energy is the average energy of an external 
electron. That is, 

f f E * <£ 0 (E\fi)dE’ dfl 

o E> 

E = - ( (50) 

/ / * o <E\?f) dE' dfi 

'r* E 1 

Q r qi 

According to Evans , the total source in MeV/g 
is approximated by 


S=HZ 0 E 2 / / * 0 <E', 0) dE' dn , 

- T7- f 


(51) 


where H is a constant and Z Q is the shield material 
atomic number. Berger and SeltzerP] tabulate H 
for a cosine law electron source as a function of 
incident electron energy and shield thickness. 
Because the variation in the table is not great, 
a typical value 4 x 10"4 can be used for this 
approximation. 

Second, assume that the source plane is 
located at half the extrapolated range at ¥. The 
bremsstrahlung dose is given by: 

D(Z)=KM E SB[E,jU m (E)Z*]E 1 tt m (E)Z*] , ( 52 ) 


For p E , the energy absorption coefficient in the 
receiver, 0.033 cm 2 /g is used, which is an upper 
bound on the coefficient for water in the range 
above 0.1 MeV. The source-receiver distance is 
given by: 


Z*=Z-R(E)/2 Z>R(E)/2 . (53 . 

Hiven 8 bv nt - ^ i UUdUp faCtor in is 

y> B(E,m (E) Z*) = 
m 

1.0+ 26.47 M m (E) Z*! 1, 161 E < 0. 1 (54) 

1.0+ 1. 827 A‘ m (E) Z* O.lsEs 0.2016 

1.0+ 1.253 V I) Z* E -0, 2354 o. 2016 < e s 2 . 0 

1-0+ 1.528 M m (E) Z*i"°- 522 I >2 . 0 
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Simplifying somewhat for an aluminum shield and a 
water receiver yields 

D(Z)=2.7xlO ,2 E 2 B[E,M p1 ( E ) z *l 

Ej [n (E) Z* 1 / J ♦, (E’.OJdB'dO (55) 

I 11 _ r ' 

n 

Figure 21 shows the result of a bremsstrahlung dose 
calculation using Equation (55); Fig. 22 shows a 
similar calculation assuming a monodirectional 
source. As can be seen in comparisons with results 
of more accurate calculations, the approximation 
yields reasonable order of magnitude estimates that 
are generally conservative. 


CONCLUSION 

Although the results of the different electron 
and bremsstrahlung dose calculations developed at 
MSFC give results in relatively good agreement with 
each other, comparison with independent calculation 
is useful in pointing up areas of possible weak- 
ness. Wayne Scott of Oak Ridge National Laboratory 
has made such calculation using a Boeing program 
called Charge and a Douglas program called BEP as 
well as an older version of the MSFC program using 
Berger’s transmission and reflection data^lO] . The 
spectrum used was 



•HIILO THWKNIM Worn* - ALUMINUM) 


rlcum: III OJWAAIliM or MEMSSTRAMUMc: DUSES CALCULATED ASSUMING 
A HALF-SPACE ISOTROFIC SOURCE AND INCIDENT SPECTRUM [N 



SHIELD THICKNESS (,'cm 2 ■ ALUMINUM) 


fICUKE 22 : COMPARISON OF BREMSSTRAHLUNG DOSE CALCULATED ASSUMING 

A MONODIRECTIONAL SOURCE DIRECTED NORMALLY INWARD AND 
A HALF-SPACE ISOTROPIC ELECTRON SPECTRUM OF THE FORM 
♦ 0 (E) - P «p(-PE)./c» 2 - M.V IN THE COMPLEX MSFC 
CALCULATION [EQUATION (47)] AND USING EQUATION (52) 
BUT WITH AN EXPONENTIAL ATTENUATION KERNEL 


-0. sfSE-0. 055E 2 , f ,(56) 

$ 0 (E) = 3. 88 e cnr-MeV-sec 


which was renormalized to a total flux of unity. 
Table 3 gives the calculated electron dose in units 
of rad/h behind an aluminum shield on a tissue 
receiver for a normally incident beam. The MSFC 
results shown are for the new program using 
Berger’s energy deposition coefficients and for the 
straight-ahead approximation. Both the MSFC pro- 
gram results are higher than the Charge and BEP 
results. The calculation using Berger’s data is 
believed to be higher because the receiver had 
infinite backing behind it, whereas the other pro- 
grams assume no backing. The straight-ahead cal- 
culation was higher because the extrapolated range 
was used in describing the transport. 

Scott also did some spectral calculation at 
0.5 and 1.0 g/cm^ for the same incident spectrum 
using Charge, BEP, and Berger's Monte Carlo pro- 
gram, ETRAN. Figures 23 and 24 show these results 
plus spectra calculated with the MSFC straight- 
ahead program; it can be seen why the straight- 
ahead program yields higher doses. For the low- 

energy end of the spectrum consisting of particles 
that have lost the largest fraction of their energy 
and have been scattered the most, the straight- 
ahead program overestimates the flux by as much as 
a factor of three compared to ETRAN. For higher 
energies there is relatively good agreement with 
ETRAN. The overestimate of the low-energy com- 
ponent of the electron flux will be reflected in 
the bremsstrahlung calculation as an overestimate 
of the low-energy component of the bremsstrahlung 
source and in turn an overestimate of the brem- 
sstrahlung dose, especially behind the thin shield 
where the low-energy bremsstrahlung has not been 
attenuated greatly. 



0 1 2 3 4 5 6? 

ELECTRON ENERGY (MiV) 

figur* xa. ■ ktctron behind u toml -inflnlt« slob uf 

iiLumiQum of thick/* «♦ ll.ft g/cm*. 
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In conclusion, the program using Berger's 
energy deposition data should be preferred for dose 
calculations, because in most practical problems 
the receiver does have effectively infinite backing 
and because the program is relatively easy to use. 
Where an electron spectrum at an internal point in 
a shield is needed, the straight-ahead approxi- 
mation can be used to obtain a conservative esti- 
mate of the differential flux. For first-order 
calculations an estimate of bremsstrahlung dose 
can be obtained by using Equation (55) where the 
incident electron spectrum falls off exponentially. 
Where an accurate calculation is needed, Equation 
(47) with a half-space isotropic attenuation kernel 
should be used. More work needs to be done here 
in determining the actual source angular distri- 
bution. 



FlKUr* 2*. Tn<r»initted olectron tpectru behind :i Mml-inlinilQ «l:tb of 
iilumioum uf thickness 
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